Electronic structure and Magnetism in BaMn 2 As 2 and BaMn 2 Sb 2 
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We study the properties of ThCr2Si2 structure BaMn2As2 and BaMn2Sb2 using density func- 
tional calculations of the electronic and magnetic as well experimental measurements on single 
crystal samples of BaMn2As2. These materials are local moment magnets with moderate band 
gap antiferromagnetic semiconducting ground states. The electronic structures show substantial 
Mn - pnictogen hybridization, which stabilizes an intermediate spin configuration for the nominally 
d 5 Mn. The results are discussed in the context of possible thermoelectric applications and the 
relationship with the corresponding iron / cobalt / nickel compounds Ba(Fe,Co,Ni)2As2. 

PACS numbers: 71.20.Lp,72.15.Jf,75.10.Lp 
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I. INTRODUCTION 

Requirements for a thermoelectric material to be used 
in waste heat recovery include a high figure of merit, 
ZT = aS 2 T '/ k, preferably unity or higher, and materials 
cost compatible with the application. Here 5* is the See- 
beck coefficient (thermopowcr), a is the electrical conduc- 
tivity, and k is the thermal conductivity. High ZT cannot 
be obtained without high thcrmopowerii Current materi- 
als for intermediate temperatures typical of requirements 
for automotive exhaust heat conversion contain substan- 
tial amounts of Te, as in e.g. PbTe. However, Te is a 
rare element whose supply is mainly as a byproduct of 
Cu production and whose cost is expected to sensitive 
to usage. In particular, this may lead to difficulties in 
applications that require large volumes of Te based ther- 
moelectric material. Therefore it is of interest to identify 
alternative materials based on common inexpensive ele- 
ments. Among metals and metalloids these include alka- 
line earths (Mg,Ca,Sr,Ba), mid-3 d transition elements, 
especially Mn and Fc and pnictogens (As,Sb). Wang and 
co-workers^ recently reported measurements of the trans- 
port properties of single crystals of BaMn2Sb2 grown 
from Sn flux. They report a room temperature ther- 
mopower of ~ 225 /iV/K, but accompanied by a low elec- 
trical conductivity leading to low ZT. The conductivity 
increases with temperature, but this is accompanied by 
a decreasing thermopowcr. This suggests intrinsic semi- 
conducting behavior, which is consistent with activated 
fits of the conductivity. Interestingly, a transition was 
found from a low activation energy of ~ 0.1 eV (T <470 
K) to a higher activation energy region (~ 1.1 eV) at 
T >650 K, possibly reflecting a cross-over from an ex- 
trinsic to an intrinsic regime. However, doping studies 
have not yet been reported. 

Recent discoveries of superconductivity^ in ThCr2Si2 
structure pnictides has led to renewed interest in this 
large class of compounds^ Among the compounds 
AM2PTI2, yl=Ca,Sr,Ba,Eu, Pn= As, high temperature 



superconductivity is found for M=Fe, and an apparently 
different superconductivity for M=Niji&£ The M=Co 
compound, BaCo2As2 is near ferromagnetism^S while 
CaCo2P2 is reported as having ferromagnetically ordered 
Co planes, which are stacked antiferromagneticallyjii In- 
terestingly, unconventional superconductivity is found in 
Ba(Fei_ a: Co a: )2As2 for a substantial range of x around 
0.1 Turning to the electronic structures, the com- 
pounds with M=Fe,Co and Ni show roughly similar 
shapes of their densities of states, with the main dif- 
ference being the position of the Fermi energy, due to 
the different electron counts. These electronic struc- 
tures show moderate hybridization between the transi- 
tion metal d bands and the As p bands, with the main 
As p bands ~ 2-3 eV below the Fermi energy, Ep and 
approximately ~ 10-20% As characater in the d bands. 
Thus the As occurs nominally as As 3 ~. Interestingly, 
even though ThCr2Si2 is an very common structure type, 
the Pn=Sb compounds, namely BaFe 2 Sb 2 , BaCo 2 Sb2 
and BaNi2Sb2 are not reported. On the other hand, the 
Mn compounds, BaMii2P2, BaMn2As2 and BaMii2Sl:>2 
are known phases as are the corresponding Ca, Sr and 
Eu materials, though with structural distortions in some 
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From a practical point of view, there are several re- 
ports of remarkably high thermopowers accompanied by 
metallic conduction in the Fe based superconducting 
phasesi 18 i 19 ' 20 i 21 ' 22 These materials have relatively nar- 
row manifolds of Fe d bands near the Fermi energy. 
These result from the combination of moderate Fe-As 
hybridization and expanded Fe-Fc distance as compared 
with bulk Fe. In any case, the values of S in the Fe com- 
pounds are too low by a factor of ~ two compared with 
practical thermoelectrics, and are high at low tempera- 
tures, inapplicable to waste heat recovery. As mentioned, 
BaMn2Sb2 is reported to be a semiconductor i2ill We note 
that BaMn2Sb2, BaMn2As2 and BaFe2As2 are composed 
of common, inexpensive elements, which as mentioned is 
an important consideration for some thermoelectric ap- 
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plications. Key questions concern the relationship be- 
tween BaMn 2 As2 and BaMn 2 Sb2, and that between the 
Mn and Fe compounds, as well as the related issue of 
whether the conductivity in the Mn compounds can be 
improved by doping, while at the same time maintaining 
high thcrmopowers. 



II. THEORETICAL METHODS 

The purpose of this paper is to further elucidate the 
electronic properties of these Mn compounds in relation 
to the Fe phases, in particular by comparing BaMn2As2 
with BaFe 2 As 2 , and with BaMn 2 Sb 2 , which was recently 
studied by Xia and co-workers.— We find that the Mn 
compounds are rather different from the Fe, Co and Ni 
compounds in that there is much stronger hybridization 
with the pnictogen p states. This leads to an electronic 
structure that cannot be regarded as derived from that 
of the Fe compound with a different electron count . Fur- 
thermore, in contrast to the metallic character and itin- 
erant magnetism in the Fe phases, we find local mo- 
ment, semiconducting behavior in both BaMn 2 Sb 2 and 
BaMn 2 As 2 . 

Electronic structure calculations were performed 
within density functional theory. We did calculations 
using both the local density approximation (LDA) and 
the generalized gradient approximation (GGA) of Perdew 
and co-workers^ and also using both the general poten- 
tial linearized augmented planewave (LAPW) method 2 ^ 
similar to our earlier calculations for the Fe-bascd 
materials ; 25 ' 26 and an ultrasoft pscudopotcntial method 
implemented in the Quantum Espresso packaged For 
the LAPW calculations, we used two different codes, 
an in-house code and the Wien2k code^ and com- 
pared results, finding no significant differences. We used 
LAPW sphere radii of 2.1 Bohr for the elements other 
than Ba in BaMn 2 As 2 and 2.2 Bohr for all elements in 
BaMn 2 Sb 2 . Two different sphere radii were used for Ba 
in BaMn 2 As 2 : 2.35 Bohr with Wicn2k and 2.2 Bohr in 
the other LAPW calculations. In all cases, converged ba- 
sis sets, corresponding to Rmmk m ax=9 with additional 
local orbitals 2 ^ were used, where R m in is the minimum 
LAPW sphere radius and k max is the planewave cut- 
off. Convergence of the Brilloiun zone sampling was 
checked by directly calculating results using different 
grids. We used the literature values of the lattice parame- 
ters, a=4.418A, c=14.200A, for BaMn 2 Sb 2 and a=4.15A, 
c=13.47A, for BaMn 2 As 2 with internal parameters, za s 
and zsb determined by LDA total energy minimization 
in the low energy antiferromagnetic state. The resulting 
values were z,sf,=0.3594 for BaMn 2 Sb 2 and z^ s =0.3524 
for BaMn 2 As 2 . The value obtained for BaMn 2 Sb 2 is 
in reasonable accord with the experimental valued of 
0.3663. 

Transport calculations were performed within the con- 
stant scattering time approximation using the BoltzTraP 
code<^ For consistency we focus here on results obtained 



within the LDA with the LAPW method, but we note 
that the conclusions are robust with very similar results 
and conclusions emerging from the GGA planewave cal- 
culations. We also note that our electronic structure for 
BaMn 2 Sb 2 is quite similar to that reported previously by 
Xia and co-workers J^- 



III. SYNTHESIS 

Experimental characterization of BaMn 2 As 2 was per- 
formed using single crystals. These were grown out of 
MnAs flux. The typical crystal sizes were ~ 5x3x0.2 
mm 3 . High purity elements (> 99.9%) were used in the 
preparation of the samples. The source was Alfa Aesar. 
First, MnAs binary was prepared by placing mixtures of 
As, Mn pieces in a silica tube. These were reacted slowly 
by heating to 300 °C (25 °C/hour, dwell 10 hours), and 
to 600 °C (15 °C/hour, dwell 75 hours). Then a ratio of 
Ba:MnAs = 1:5 was heated in an alumina crucible for 15 
hours at 1230 ° C under partial atm argon. This reaction 
was cooled at the rate of 2 °/hour, followed by decanting 
of the MnAs flux at 1120 °C. 



IV. CRYSTAL STRUCTURE, TRANSPORT 
AND SPECIFIC HEAT 

The crystals were malleable but well-formed plates 
with the [001] direction perpendicular to the large faces. 
The structural identification was made via powder x-ray 
diffraction using a Scintag XDS 2000 6-0 diffractomcter 
(Cu K a radiation). The cell parameters were refined us- 
ing least squres fitting of the measured peak positions in 
the range 29 = 10 - 90° using the Jade 6.1 MDI pack- 
age. The resulting lattice parameters of BaMn 2 As 2 were 
a=4.1633(9) A, c=13.448(4) A, in the ThCr 2 Si 2 struc- 
ture (14/mmm, Z~2). The cell volume for BaMn 2 As 2 is 
then 233.087(4) A 3 , which is in accord with prior reports 
and much larger than BaFe 2 As 2 (204.567(2) A 3 ) (Rcf. 
GJ) or BaCo 2 As 2 (197.784(1) A 3 ) (Ref. [13). 

Temperature dependent electrical resistivity measure- 
ments were performed on a Quantum Design Physical 
Property Measurement System (PPMS) , measured in the 
afr-plane above 5 K. For BaMn 2 As 2 , the resistivity, p de- 
creases with increasing temperature up to ~ 150 K, above 
which it gives a decreasing, metallic-like dependence (Fig. 
[1}. The resistivity at room temperature is 165 mfi cm. 

Specific heat data, C P (T), were also obtained using the 
PPMS via the relaxation method. Fig. [3] gives the tem- 
perature dependence of the specific heat. There are no 
features suggesting a phase transition in the specific heat 
up to 200 K. Below ~6K, C/T has a linear T 2 depen- 
dence (inset of Fig. [5]). The fitted Sommerfeld coefficient, 
7=0.79(2) mJ/(K 2 mol), or ~ 0.4 mJ/(K 2 mol Mn) is 
near zero. For comparison the value fo BaFe 2 As 2 , which 
has a spin density wave instability gapping most of the 
Fermi surface is 3.0 mJ/(K 2 mol Fe)^ and the value for 
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FIG. 1: Temperature dependence of the in-plane resistivity 
of BaMn 2 As 2 betrween 1.8 K and 380 K. 



BaCo 2 As 2 is 20.8 mJ/(K 2 mol Co).^. In fact, it is likely 
that the small observed 7 arises from a small amount of 
MnAs flux included in the sample and is therefore not 
intrinsic to BaMn2As2- This is based on squid measure- 
ments, where we find a small signal at 310 K, correspond- 
ing to the ordering temperature of MnAs. The value for 
the Debye temperature is Qd ~ 280 K, estimated above 
150 K, using the calculated values of the T/Op depen- 
dence of the Debye specific heat model. Thus our trans- 
port and specific heat measurements are consistent with 
a small band gap semiconductor. They are not consistent 
with either a large gap insulator (e.g. a Mott insulator) 
or a metal with significant carrier density. 



V. ELECTRONIC STRUCTURE AND 
MAGNETISM 

We begin the discussion of the density functional in- 
vestigation with a summary of our calculated results 
for BaMii2Sb2, which as mentioned are similar to those 
obtained previously by Xia and co-workers ^ although 
in general the LAPW method is more accurate than 
the linear muffin tin orbital method employed in that 
work. Calculations were performed for a checkerboard 
antifcrromagnctic ordering (AF-C), which is the lowest 
energy state, as well as a ferromagnetic ordering (F), 
and zone boundary antifcrromagnctic orders alternating 
chains of like spin Mn (as in the spin density wave state 
of BaFe2As2), and a non-spin-polarized (NSP) case. The 
two cells with chain ordering are denoted AF-S1, which 
consists of chains stacked antifcrromagnctically along the 
c-axis direction, and AF-S2, with ferromagnetic stack- 
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FIG. 2: Temperature dependence of the specific heat for 
BaMn2As2. The inset shows C/T vs. T 2 and a linear fit 
between 1.8 K and 3.5 K 



ing. The energetics are summarized along with those of 
BaMn2As2 in Table HI We find a strong tendency to- 
wards moment formation, which results from the strong 
Hund's rule coupling on the nominal Mn 2+ ions in the 
compound, with stable moments independent of the mag- 
netic order. Integrating within the Mn LAPW sphere 
(radius 2.2 Bohr) we obtain moments of 3.47 hb for the 
ferromagnetic (F) ordering and 3.55 (1b for the G-type 
checkerboard AF-C ordering, which is the lowest energy 
state (n.b. G-type ordering is the experimental ground 
state at least for the phosphide) pi£ Including contribu- 
tions from all atoms in the unit cell, the F ordered spin 
magnetization is 3.76 [13 on a per Mn basis, reflecting 
a small parallel contribution from the Sb. Significantly, 
this is considerably reduced from the nominal 5 /jlm ex- 
pected for high spin Mn 2+ . This reduction reflects strong 
hybridization between Mn d and Sb p states in this com- 
pound. In general, very strong hybridization is needed to 
effectively compete against the very strong Hund's inter- 
action in the half-filled d shell of Mn 2+ though we note 
that strong hybridization was previously identified in the 
phosphides by Hoffmann and Zheng and by Gustenau 
and co-workersi 31 i 32 ' 33 Also, associated with the strong 
hybridization we find a very large energy difference be- 
tween the ferromagnetic and antifcrromagnctic ordcrings. 
This amounts to more than 0.29 eV/Mn in the LDA and 
implies a very high magnetic ordering temperature well 
above room temperature. This would be similar to the 
phosphide, which is reported to have an ordering temper- 
ature above 750KJ£ For comparison the corresponding 
energy for bec Fc (Curie temperature, 7c=1043K) is 0.4 
eV, with opposite (ferromagnetic) sign^i 

In fact, strong, spin dependent hybridization is seen 
in the electronic structure, represented by the electronic 
density of states and projections (Fig. [3]). As may be 
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TABLE I: Magnetic energies in eV of BaMn 2 Sb 2 and 
BaMn2As2 as obtained within the LDA on a per formula unit 
(two Mn) basis. The non-spin-polarized case, denoted NSP, is 
set as zero. F denotes ferromagnetic order, while AF-C is the 
checkerboard G-type antiferromagnetic ordering, and AF-S1 
and AF-S2 are stripe- like orders (see text). 
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FIG. 3: (color online) Electronic DOS and projections onto 
the LAPW spheres for antiferromagnetic BaMn 2 Sb2. The 
plot shows the total DOS above and below the axis, the Mn 
majority spin projection above and the Mn minority spin pro- 
jection below. The remaining (non-Mn) contribution to the 
valence bands is mainly from Sb p states. 



noted, the antiferromagnetic structure is found to be a 
small band gap semiconductor, E g ~0.21 eV in the LDA. 
The majority spin Mn d states overlap in energy with 
the Sb p states and are strongly mixed with them. The 
minority spin Mn d states are above the main Sb p DOS, 
and are therefore less strongly mixed with the Sb states. 
Nonetheless, there is sufficient minority spin Mn d - Sb p 
hybridization to lead to noticeable Mn d character in the 
occupied valence band states. This is the main source of 
the moment reduction, corresponding to a Mn d electron 
count closer to Mn 1+ than Mn 2+ . Finally we note that 
such strong spin dependent hybridization while leading to 
strong superexchange interactions is highly unfavorable 
for carrier mobility in an antiferromagnet,— since at low 
T hopping between nearest neighbor sites with opposite 
spin will be suppressed and at higher T spin disorder will 
lead to strong scattering. 

We now discuss our experimental and calculated re- 
sults for BaMn2As2. The calculated magnetic energies 
(Table [l| show strong local moment magnetism with a 
large exchange, similar to BaMn 2 Sb 2 , with the exception 
that the energies of the magnetically ordered states are 



not as low compared with the non-spin-polarizcd case, as 
those in the antimonide. This reflects lower moments due 
to even stronger hybridization. The LDA total spin mag- 
netization for the F ordering is 2.9 /xs/Mn, i.e. ~ 0.9 [Ib 
lower than in the antimonide. The moment inside the 2.1 
Bohr LAPW sphere is 3.20 hb for the AF ordering and 
2.74 /j,b for the F ordering. Nonetheless, the ordering 
energy remains very high, again reflecting strong spin- 
dependent hybridization. The energy difference between 
the AF and F orderings is 0.33 eV / Mn in the LDA. 
The corresponding GGA value of 0.38 eV is similar, the 
main difference between the LDA and GGA results be- 
ing that the moment formation is stronger in the GGA, 
leading to moments in the LAPW spheres of 3.07 [ib for 
F ordering and 3.35 [ib for AF ordering. Turning to the 
AF-S1 and AF-S2 orders, the energies are slightly lower 
than the average of the F and AF-C states. A simple 
nearest and next nearest neighbor superexchange picture 
is probably an oversimplification for these narrow gap 
semiconductors. However, the results indicate that the 
nearest neighbor interaction is dominant over further in- 
teractions and that the c-axis interaction is considerably 
weaker than the in-plane interaction and that it is anti- 
ferromagnetic. 

The calculated LDA electronic density of states of 
BaMn 2 As2 is shown in Fig. [4] for both F and AF order- 
ings. The low energy AF state again has a small semi- 
conducting gap of 0.1 eV in the LDA and 0.2 eV in our 
GGA calculations. These low values are consistent with 
the activation energies from resistivity data below 470 K 
obtained by Wang and co-workers^ Thus we interpret 
their low temperature value of ~ 0.1 eV as the intrinsic 
bulk gap for the AF ordered state. Within this scenario, 
the larger activation gap at high T would be associated 
with Anderson localization arising from magnetic disor- 
der above the AF ordering temperature. 

Strong spin-dependent hybridization is seen, and in 
particular there is even more minority spin Mn charac- 
ter mixed with the pnictogen derived valence bands than 
in the antimonide. This stronger hybridization is consis- 
tent with the general trend that As has a greater ten- 
dency towards bond formation with transition elements 
than Sb and is reflected in the higher energy difference 
between the F and AF magnetic orderings in BaMn2As2 
than in BaMn 2 Sb2. We emphasize that this behavior is 
very different from that found in BaFe2As2, BaCo2As2 
and BaNi2As2, which have some d - p hybridization, as in 
an oxide, but are much more ionic3 ' 10 i 26 The distinct Mn 
- As bonding in this structure also explains why Mn is 
not an effective dopant for the Fe-based superconducting 
compounds, while Co and Ni are^ 



VI. THERMOPOWER 

As noted, antifcrromagnctism with strong spin de- 
pendent hybridization is not favorable for high mobility 
conduction. Nonetheless, the chemical flexibility of the 
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FIG. 4: (color online) Electronic DOS and projections onto 
the LAPW spheres for ferromagnetic ordered (top) and anti- 
ferromagnetic BaMn2As2 (bottom). The plots show majority 
spin above the axis and minority spin below. The remaining 
(non-Mn) contribution to the valence bands is mainly from 
As p states. 
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FIG. 5: (color online) LDA constant scattering time approx- 
imation Seebeck coefficient in the basal plane of AF ordered 
BaMn2As2. The doping levels are in carriers per formula unit 
(2 Mn). 
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FIG. 6: (color online) LDA constant scattering time approx- 
imation c-axis Seebeck coefficient of AF ordered BaMn2As2. 
The doping levels are in carriers per formula unit. 



ThCr 2 Si2 structure suggests that high doping levels may 
be achievable. In that case, if the thermopower remains 
high, as in the oxide thermoelectric Na^OoO^, high ZT 
may result in spite of low mobility^ Therefore we calcu- 
lated the Seebeck coefficients from the LDA band struc- 
ture as a function of doping level and temperature in the 
constant scattering time approximation. The results are 
shown in Figs. [5]and[5J The anisotropics are different for 
the valence and conduction bands. In the basal plane, the 
hole doping yields higher S than electron doping, while 
an opposite trend is found in the conduction bands. This 
difference is due to different band anisotropics. The con- 
stant scattering time transport anisotropy a zz ja xx esti- 
mated from the band anisotropy is 1.9 for the valence 
bands (p-type doping) and 0.5 for the conduction bands 
(n type). These modest anisotropies mean that the mate- 
rial is effectively a three dimensional intermetallic com- 
pound rather than a quasi-2D system. Turning to the 
values of S we observe that high values compatible with 
high ZT are only found for modest doping levels and not 
for high metallic doping levels such as those in Na^OoO^ 
(0.01 carriers per formula unit is 8.6xl0 19 cm~ 3 ; the car- 
rier concentration in thermoelectric Na 2 .Co02 is more 
than one order of magnitude higher) Considering the 
likelihood of carrier localization in this magnetic mate- 
rial, this would initially seem not favorable for applica- 
tion of BaMii2As2 as a thermoelectric for waste heat re- 
covery. 



VII. SUMMARY AND CONCLUSIONS 

To summarize, we characterize BaMn2As2 and 
BaMn2Sb2 as small band-gap, local moment antiferro- 
magnetic semiconductors. We find a very strong spin 
dependent hyridization. This leads to high exchange en- 
ergies, which are consistent with high magnetic ordering 
temperatures. The lowest energy state found is the G- 
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type checkerboard antiferromagnetic state for both com- 
pounds. We note that the bonding, electronic structure 
and magnetic properties arc very distinct from those of 
the corresponding Fe, Co and Ni compounds, which may 
explain why Mn leads to carrier localization rather than 
effective doping in BaFe2As2 and related materials. 

The strong spin dependent hybridization is not favor- 
able for high carrier mobility. Nonetheless, the high thcr- 
mopowers do suggest that experimental doping studies 
should be performed to determine the maximum ZT in 
this material. 
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